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■ Differential spectra of coherent photon bremsstrahlung in relativistic 

heavy ion collisions are calculated within various schematic models of the 
projectile-target stopping. Two versions of the degradation length model, 
based on a phenomenological deceleration law, are considered. The simple 
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£T) ■ shock wave model is studied analytically. The predictions of these models 
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agree in the soft photon limit, where the spectrum is determined only by the 
final velocity distribution of charged particles. The results of these models 
in the case of central Au+Au collisions at various bombarding energies are 



compared with the predictions of the microscopic transport model UrQMD. 
It is shown that at the AGS energy the coherent photon bremsstrahlung 

tx! - 

j_j ■ exceeds the photon yield from 7r°-decays at photon energies to < 50 MeV. 
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I. INTRODUCTION 

Real photon bremsstrahlung emitted in the course of a heavy ion collision may provide 
information about the underlying dynamics and the stopping law of nuclear matter. To 
understand the sensitivity of photon spectra to the dynamics of a heavy ion collision, it is 
instructive to compare the predictions of models based on different dynamical scenarios 
of the collision process. In earlier publications the coherent bremsstrahlung was studied 
in phenomenological |[]-f| as well as microscopical || models. In the present paper we 
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concentrate on the comparative analysis of different models, which has not been done in 
the above mentioned works. 

In section [TT] some general properties of bremsstrahlung spectra common to all models 
are obtained from basic considerations about the nature of coherent bremsstrahlung. 
section JTJ introduces different deceleration models: a shock wave model and a degradation 
length model both for point-like and extended nuclei. In addition we employ a microscopic 
transport model which is not restricted to the initial deceleration stage. The comparison 



of these models is performed in section [TV. 

A serious problem in observing the coherent photon bremsstrahlung is the large back- 
ground of 7T° — > 77 decays. It is especially strong at ultrarelativistic bombarding energies 
(AGS and higher). To overcome this problem it was suggested to study experimen- 
tally sufficiently soft photons (with energies less than several MeV in the case of RHIC) . 
According to our calculation the contribution of the ir° decays is indeed relatively small 
at low photon energies. This problem will be studied in more detail when we compare 
the models in section |TV| . Conclusions and discussion will be given in section M. 



II. GENERAL DISCUSSION 

In this section some general model-independent properties of bremsstrahlung spec- 
tra are discussed. According to the Larmor formula the energy spectrum of the 
bremsstrahlung is proportional to the squared absolute value of the Fourier transformed 
force f(t) acting on a charge ||10|| , dl/dw ~ \ f{uj)\ 2 . From general properties of the Fourier 
transformation one is lead to a number of conclusions: 

i) At large photon energies the radiated energy tends to zero if the force is absolutely 
integrable. In contrast, instantaneous interactions (5 forces) generate a flat energy spec- 
trum. The latter approximation is valid in the so-called soft-photon limit (u — > 0), which 
in turn can be interpreted as a limit r — > where r is the interaction time. 

2 



ii) From the scaling behavior of Fourier transforms it may be shown that the 
bremsstrahlung drops the faster in uj the longer the force acts on the charge (this is 
the analogue of the Heisenberg's uncertainty principle). As a result, the energy spectra 
of bremsstrahlung photons become harder in the collision of smaller nuclei or at higher 
bombarding energies. 

hi) If the force f(t) is essentially nonzero in a finite time interval \t\ < a, which is a 
common assumption in all models describing heavy-ion collisions far above the Coulomb 
barrier, the spectra of bremsstrahlung photons exhibit oscillations in uo. Indeed, any con- 
tinuous function on a finite interval can be approximated with arbitrary accuracy by a 
polynomial of a certain order n, so that it is sufficient to consider only forces of polynomial 
type. By subsequent integration by parts, it can be shown that the Fourier transform 
f(u) may be written as f(u) = Ae iuja + Be"™, where A, B = =f £ i3+1 ^+i ±a) and A(B) 
corresponds to the upper (lower) sign. Therefore the bremsstrahlung spectrum dl/doo 
consists of a polynomial part superimposed on an oscillating part. In this simple example 
the oscillating part has the frequency ir/a. In real physical processes the frequency of the 
oscillations depends upon details of the acting force. 

In a semiclassical approximation the radiated energy per invariant phase-space element 
is given by |T| 

d 3 I 

2(2tt) 3 — = | j;(k)f(k) | . (1) 

where / = uN, k = (u, k) is the photon 4-momentum and j^(k) is the Fourier transform 
of the classical 4-current. By using charge conservation, k lJL j IJ, (k) = 0, and neglecting 
the transverse components of j(k), which is a good approximation at high bombarding 
energies, one has 

\mnk)\ = (^j \ 3z (k)\ 2 - (2) 

From Eqs. (|1]) and (fj) it follows that 
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Let us assume that photons are produced in microscopic binary collisions of point-like 
particles which instantaneously change their velocities in each collision vertex. Then the 
explicit expression for j' 1 (k) may be written as 

f{k) = i EE * ( r^- - r^-) ^ > ( 4 ) 

where indices i and j count vertices and particles, respectively. The 4-vector is the 
space-time position of the z-th collision vertex of the j-th particle, Cj is the charge of 
the latter. It is assumed that at x — the j-th particle 4-momentum jumps from Pi-ij 
to pij . Eq. (|4]) can easily be applied to inelastic particle collisions if one extends the 
definition of the charge ej — > and allow the charge to change at a vertex. For example, 
if some charged particle is produced at the vertex i then e^ij = and 7^ 0. 

Assuming for simplicity that the particles move along the beam direction z one immedi- 
ately finds that the emission of soft photons is suppressed if the particles are reaccelerated 
in the course of the reaction. The same phenomenon takes place if charges are produced 
and reabsorbed during a heavy-ion collision. By introducing a time cut one can decom- 
pose the current into parts having fixed but opposite signs of the acceleration or into parts 
corresponding to increasing and decreasing number of charged particles, respectively. It 
can be shown that the interference of the radiation from these two current components is 
destructive for soft photons. 

The soft photon limit of (0) is independent of the collision dynamics, but depends 
only on the charged particles' velocities in the initial and final states of the reaction. In 
this case the sums over the in- and outgoing particles can be replaced by integrals over 
the initial and final velocity distributions. Let us consider the soft-photon limit of ([[]) for 
the case of symmetric central heavy-ion collisions, assuming that the charges move along 
the beam axis. For sharp initial and final velocity distributions one finds 



dl _ AZ 2 a{(5 2 - (3}) 2 sm 2 d cos 2 d 
dudn ~ (2tt) 2 (1 - 0? cos 2 tf) 2 (l - fa] cos 2 $) 2 ' !ii 

Here fa is the initial (final) velocity of charged particles, d is the cm. emission 
angle, Z is the total charge of each nucleus and a = e 2 / (Anhc) is the QED fine structure 
constant. The l.h.s. of Eq. (f|) exhibits the typical behavior of quadrupole radiation, 
which is a direct consequence of the symmetry of the considered system. 

Using Eq. (M) it can be shown that the angle of maximum radiation equals 



4ax = arccos 



U-L- — U 



with 



u± 



-A + 2(3 2 + 2(3 2 -3(3 2 (3 2 \ 3 (fi + 0} + &0}^ 



Clearly for fa >0.9 ?9 max depends quite weakly upon j3f (if /?/ is not too close to fa). To a 
good accuracy it may be assumed that (3f ~ 0, which leads to 

$max = arcctg7i (6) 

with 7j = (1 — (3 2 ) 1 / 2 . Substituting -d = $ max into fl5|) yields the following estimate for the 
height of the maximum in the photon angular distribution 

dl 



lim 



o duodVt 



Z 2 a(P 2 -P 2 ) 2 (2-p 2 ) 2 
x (27r) 2 (l-/3 2 )(2-/? 2 -/3 2 ) 2 U 

If (3f ~ fa the exact value of $ max becomes smaller. In the case /3/ ~ fa ~ 1 one finds 
"^max ~ l/(v37i)- At large bombarding energies the photons emitted in a heavy-ion 
collision are strongly peaked in forward and backward directions and therefore practically 
do not overlap in phase space. In this case the interference term between the projectile and 
target photons becomes negligible and the quadrupole spectrum may be well approximated 
by the sum of the dipole spectra generated by each nucleus. For pure dipole radiation we 
find $ max = arccos fa for Pf — and i? max ~ l/(\/37i) for fa ~ /3j ~ 1. Striclty speaking, 



however, using the approximation of dipole radiation in (symmetric) heavy-ion collisions 
is not correct. 

The same conclusion about the low sensitivity of the soft photon angular distribution 
also holds for broad final velocity distributions, provided that the mean final velocity 
is smaller than fy. In this case only the amplitude of the radiation is sensitive to the 
degree of stopping. This discussion shows that for a large range of nuclear stopping soft 
bremsstrahlung alone can not serve as a good tool to measure the impact parameter f2j. 
In the soft photon limit only the amplitude of the bremsstrahlung spectrum is affected by 
both the involved charge and the final velocity distribution. The investigation of the soft 
photon spectrum (Eqs. (|5|) and ([F])), however, shows that with increasing bombarding 
energy the sensitivity of the amplitude on the final velocity becomes extremely weak for 
a large range of final velocities. For example, in the case of a Au+Au collision at RHIC 
energies (yi = artanh/3j ~ 5.4) the amplitudes for the final cm. rapidities yf — and 
Uf = 3 differ by less than 5%. 

As is well known, scattering of particles exclusively along the beam direction excites 
only the m = multipole components of the radiation field. This implies that no photons 
are emitted in and transverse to the beam direction. Of course, transverse scattering or 



the loss of the projectile-target symmetry of the current (e.g. in peripheral collisions |T2 



excite components with m^O. As a consequence, nonvanishing radiation at $ ~ and 
7r/2 will appear ||. 

Let us now consider the influence of the nuclear charge form factor on the photon 
spectra. At large bombarding energies particles move practically along the beam axis, and 
transverse sizes of the colliding nuclei remain almost unchanged during the most violent 
initial stage of a heavy-ion collision. It is assumed for simplicity that the transverse form 
factor F(k±) can be factored, yielding approximately 



where (dl /dujdVL)^ denotes the spectrum obtained for slab-like nuclei. This formula is well 
justified only for a certain class of charge density distributions (see e.g. Sect. |III A|) . It 
becomes exact at very large bombarding energies, where F(k±) is given by the transverse 
Fourier transform of the nuclear thickness function M. The presence of the nuclear form 
factor F(k±) suppresses the photon yield at transverse photon momenta k±>l/R, where 
R is the nuclear radius. At fixed photon energy u, this suppression becomes stronger 
for larger angles of emission with respect to the beam axis. As a result, with increasing 
bombarding energy structures of dl/dudQ are shifted to smaller emission angles. 

Another feature common to all models is that the compression of nuclear matter at 
intermediate stages of a heavy-ion reaction is correlated with the characteristic time of the 
collision. In general, smaller collision times correspond to larger intermediate densities 
and to smaller slopes of the spectra with respect to uj. Therefore, measurement of these 
slopes may provide some information about the maximal compression of matter in a 
heavy-ion collision. 



III. DYNAMICAL MODELS OF NUCLEAR COLLISIONS 

A. Shock wave model 

In the first model it is assumed that two plane shock waves moving with constant cm. 
velocities, ±/? S h, develop at the moment of the first contact of the projectile and target 
nuclei. Within the shock wave model (SWM) it is postulated that the cm. velocity of 
nuclear matter is zero behind the shock fronts, i.e. at \z\ < (3 s ^t. It is assumed that the 
initial nuclei are cylinders with effective radius R and length 2R/ r y i . One can write the 
relation QQ^i—TtR 2 = A, where q$ — 0.17 fm -3 is the normal nuclear density and A is the 
atomic number of the colliding nuclei. In the case of Au nuclei, R ~ 0.98A 1 / 3 fm ~ 5.7 fm. 
This is slightly less than the geometrical value R g ~ 1.12A 1 / 3 fm ~ 6.5 fm. 
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Before the collision, at t < 0, the longitudinal component of the electromagnetic 
current can be written as 

Jz(t,r) = e -jQ 0li PiQ{R - r ± ) [e(j3it - z) Q(z - z x {t)) (9) 

Here 9(x) = |(1 + sgnx) , and Z\{t) = —2R/ji + Pit, r±(z) are the transverse (longi- 
tudinal) coordinates of r with respect to the beam axis. The first and the second term 
in the square brackets of Eq. (H) correspond to the projectile and target nucleus, respec- 
tively. The shock fronts, appearing at t = 0, reach the rear sides of the colliding nuclei 
at t = r = Tfjj^p h ) ■ At t > t compressed matter starts to expand into the vacuum. We 
neglect the photon production at this less violent stage of the reaction, assuming that 
j z (t,r) vanishes fot t > r. The explicit expression for j z at the intermediate stage, i.e. 
for < t < t, is given by Eq. ([]) with the replacement Q(Pit — z) — > Q(z + P s ht) . 

In accordance with Eq. ([!]), the spectrum of bremsstrahlung photons is determined by 
the Fourier transform of the 4-current j^{k) . After introducing the transverse density 
form factor of the initial nuclei 

one arrives at the expression 

I — e i^-r(l+/3 sh cos0) e iuiT (l-/3 sh cos 0) _ ^ 



3z (k) = ZeF(k ± )Pi 



'ID 



(UJ — Pi k z ) UJT (1 + P sh COS 9) (u + Pi k z ) LOT (1 — /9 s h cos^ 

Equivalently, the Fourier transformed current can be calculated directly using Eq. 
assuming homogeneously charged nuclear matter and replacing the sum by an integral. 
In this calculation the space-time positions of the vertices Xij coincide with the positions 
of the shock fronts. 

Using Eqs. (^), (|TTJ) yields the final result (the respective formula given in Ref. |I3| is 
incorrect) 
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dJ ( Z Pi sin 6 \ 2 2 . f sin 2 ^r(l + /3 sh cosg)/2] 

dudtt a { 71L0T j Sm ' ^ COS 6) 2 {l+(3 sh COS Of 

sin 2 [ur(l — P s h cos6 l )/2] cos (u;T/? S h cos#) [cos cur — cos (co>r/3 s h cos 6)] 
+ (l + Acos#) 2 (l-/? sh cos#) 2 + (1 -& 2 cos 2 fl)(l -/3 s 2 h cos 2 ^) 

This formula contains only one model parameter - the shock wave velocity /5 s h . It de- 
termines the duration of the collision and the achieved compression of nuclear matter. 
Strong shocks whith large energy densities of the stopped matter correspond to /5 s h <C 1. 
On the other hand, from the baryon current conservation, it can be shown that in the limit 
/3 s h — > 1 the compression of matter behind the shock front q/qo approaches its minimal 
value \J (1 + Pi) /(I — Pi). The soft photon limit clearly corresponds to lot <C 1 . It may 
be verified that the soft photon limit is reached for oj or r — > 0. In the lowest order the 
photon spectrum given by Eq. (H) is obtained, where (dl /dudfl) is equal to the r.h.s. of 
Eq. (ID with Pf — 0. As expected, in this limit the photon spectrum does not contain 
any information about the parameters of the shock waves. 



B. Degradation length model 

In the second model we treat each nucleus as a set of n test particles with the charge 
e(Z/n) moving along the beam axis. In the ground state of a nucleus the test particles are 
homogeneously distributed. It is assumed that these particles collide when their trajecto- 
ries intersect in the t, z plane. The stopping law for the colliding matter is determined by 
the rule how the slopes of the trajectories change at the collision vertices. In the model 
we have applied a simple stopping law [^4 



dz~ A' [L6) 

where z is the distance in nuclear matter covered by the test particles and A is the so 
called degradation length. 
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To implement (|T3|) into the degradation length model (DLM1) we postulate that the 
change of momentum of each test particle is determined by its number of collisions, i.e. by 
the number of trajectories crossed by the test particle during its path in nuclear matter. 
The computational procedure is based on the finite difference equation 

Az 

p{z + Az)=p{z)-p{z) — (14) 



The l.h.s. of Eq. (H) gives the momentum after the test particle collision, Az is a 
fixed value determined by the distance travelled by the test particles in the ground state 
of the nucleus. The vertex positions and velocities entering Eq. (||) are provided by the 
construction rule for the trajectories. In the limit n^oowe again treat a homogeneously 
charged matter and unphysical oscillations in the photon spectra caused by a finite number 
of test particles are washed out. The bremsstrahlung spectrum is obtained by substituting 
(D into ©. 

Figure [I] shows the trajectory representations of the SWM and the DLM1. Dynamical 
differences between these two stopping scenarios are clearly visible. One can see that 
in the SWM the stopping time is determined by /3 s h and independent of the degree of 
stopping whereas in the DLM1 this time depends on the degree of stopping. 

In the extended degradation length model (DLM2) we generalize the DLM1 and con- 
sider the colliding nuclei as homogeneously charged spheres represented by test particles. 
Binary collisions of test particles are treated in the same way as above. The final velocity 
of each test particle is again determined by the distance travelled in nuclear matter, but 
this distance now depends on the transverse separation from the nuclear center. Therefore 
spherical shapes of the colliding nuclei give rise to a broadening of the final velocity distri- 
bution and to a more complex structure of the current. On the contrary, the assumption 
of a cylindrical nuclear shape adopted by the SWM leads to a sharp final velocity distri- 
bution (Pf = 0). As discussed in Sect. II, the broadening of the final velocity distribution 
leads to almost the same shapes of soft photon spectra, but reduces the amplitude of the 
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photon yield. 

The influence of both the degree of transparency and the broadening of the final veloc- 
ity distribution due to the spherical geometry of the nuclei is demonstrated in Figure 0. 
Since stopping is effectively reduced due to the spherical shape of the nuclei the radiated 
energy is also reduced. The comparison of spectra predicted by DLM1 and DLM2 indeed 
confirms this conclusion. In accordance with the previous discussion the sensitivity of the 
photon spectra to the final velocity distribution is low for small final velocities, i.e. in the 
limit A — > 0. In the case of weak stopping (large A) the maxima of the angular distribu- 
tion shift to smaller angles. Calculations show that for A < 5 fm the angles of maximum 
radiation agree well with Eq. Within the DLM1 the heights of the maxima are well 
explained by Eq. ([F[) with the final velocity obtained from (0). Alternatively Eq. (^) may 
be used to estimate the mean final velocity in the DLM2. 

C. Transport model 

In the present paper we also consider the microscopic transport model UrQMD [ |l5|j . 
In this approach the whole history of the heavy ion collision enters the photon spectrum. 
This allows to take into account the contributions from the late stages of the collision. 

To demonstrate how the microscopic currents are treated in the model it is instructive 
to rewrite Eq. (HI) as 



where the notation 



6j if j corresponds to an incoming particle 
ij = { ' ' (16) 

-Sj if j corresponds to an outgoing particle 

is introduced. The sum runs over all charged particles j and their collision vertices i. 
Such a procedure is applied to any kind of vertex (elastic or inelastic collisions with arbi- 
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trary number of produced particles) as well as for decays of unstable particles or strings. 
Influence of mean-field potentials producing smooth bending of particle trajectories is 
neglected. Here it is assumed that hard binary collisions results in much stronger decel- 
eration of the particles. 

An important background contribution to photon spectra comes from 7r° — > 77 decays. 
In the UrQMD model the contribution of these decays can be calculated directly. In this 
calculation it was assumed that 7r° mesons decay isotropically in their rest frame. 

To reduce the statistical fluctuations of the calculated spectrum the latter was averaged 
over a large number of events. 



IV. COMPARISON OF PHOTON SPECTRA 

In this section we compare the model predictions in the case of central Au+Au col- 
lisions at various bombarding energies from GSI to RHIC. Since photon radiation in a 
symmetric heavy-ion collision is symmetric in the cm. frame with respect to $ = tt/2, 
the results are shown only for < $ < tt/2. 

Figure [5] shows the angular distribution of bremsstrahlung photons at the AGS bom- 
barding energy 10.6 AGeV. Four different photon energies are considered. The parameters 
of the models are y s h = artanh/3 s h = 1 (SWM) and A = 2fm (DLM). In the considered 
reaction such a degradation length corresponds to rather strong stopping. 

At u = 5 MeV the SWM and the DLM1 overestimate the bremsstrahlung yield as 
compared to the UrQMD calculation. The inclusion of a realistic nuclear shape, however, 
results in a good agreement of the DLM2 with the UrQMD model. This agreement favours 
a strong stopping at intermediate stages of a central Au+Au collision at AGS energy. Note 
that the degradation length A entering the model calculation A = 2 fm corresponds to 
A' ~ 5.2 fm in the laboratory frame. The latter is close to the value obtained by fitting 
pA data ]H] . 
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All simple models incorporating only the stopping phase of a heavy ion collision un- 
derestimate the radiation in forward (i? ~ 0) and transverse (i9 ~ 7r /2) direction to the 
beam. This is especially visible at larger uj's. As will be shown below, the additional 
contributions originate mainly from the late stages of the heavy-ion collision. Indeed, 
contributions from the expansion phase of the collision, and bremsstrahlung produced in 
rescatterings with newly produced particles are disregarded in the SWM and DLM. It 
should be noted that at the considered energy the transverse scattering of particles (see 
Ref. ||) gives only small contributions and cannot explain the above-mentioned discrep- 
ancy. 

We find that up to w = 50 MeV all models predict roughly the same angular distri- 
butions. In this case photon spectra are dominated by the quadrupolar radiation pro- 
duced in the course of mutual deceleration of the colliding nuclei. As discussed above, 
bremsstrahlung photons with energies less than or of the order 1/r, where r is the collision 
time, cannot resolve the detailed structure of the current. On the other hand, the simple 
models predict a more rapid decrease with u than the UrQMD model. 

Within the SWM and the DLM2 the decrease of the photon yield at uj ^> 1/r is 
determined mainly by the nuclear form factor. The shift of maximum radiation to smaller 
angles (see Sect. |T|) is clearly visible already at uj = 150 MeV. As compared to the DLM1, 
the inclusion of the form factor in the DLM2 leads to the appearance of additional maxima. 
Note that predictions of the SWM and the DLM2 are rather close, although these models 
are based on very different dynamical scenarios. The SWM predicts additional maxima 
with comparable heights and locations. It is interesting that these maxima appear in this 
model for any /3 s h ^ even for slab-like nuclei. 

At uj = 250 MeV angular distributions calculated within the simple stopping models 
and the UrQMD model differ considerably both in shapes and amplitudes of the photon 
spectra. At this energy additional maxima appear also in the UrQMD model. 
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As noted above, the spectra predicted by the SWM and the DLM decrease faster with 
respect to uj as compared to the spectrum calculated within the UrQMD model. Special 
analysis shows that such a behaviour does not originate solely from different collision 
times included in these models. Indeed, at the AGS energy the spectrum predicted by 
the SWM decreases faster with uj than the DLM2. On the other hand, the characteristic 
collision time of the DLM2 (r ~ 8fm/c) is more than twice as large as compared to the 
SWM (r ~ 3.5fm/c). A rough estimate shows that the collision time of the UrQMD 
model is of the order of 6fm/c. 

The freedom in the choice of model parameters is limited by the assumed stopping 
law. For instance, it is not possible to improve the agreement between the DLM2 and 
the UrQMD model at higher photon energies (by chosing a smaller A) without losing 
agreement for soft photons. On the other hand, in addition to the nuclear form factor, 
only the magnitude of the collision time determines the slopes of the photon spectra 
within the SWM and DLM2. 

To clarify the situation, a more detailed study of the photon spectra has been per- 
formed within the UrQMD model. Figure |4] shows the different components of the 
bremsstrahlung spectra generated in UrQMD events. A time cut at t cut ~ 6 fm has 
been introduced to separate the contributions from the early stopping phase and the late 
stages of the collision. This time cut is close to the passage time of the colliding nu- 
clei t ~ 2R/{ r )i(3i) where R ~ 7 fm is the geometric radius of the gold nucleus. The 
contribution of the stopping phase (t < t cnt ) has the familiar quadrupolar shape and 
is well understood within the SWM and DLM. Bremsstrahlung from the late stages is 
nearly isotropic and is smaller for soft photons. This is understandable since the charge- 
current varies smoothly with time at these stages of the collision. The characteristics of 
the bremsstrahlung component from t > t cut do not contradict the conclusion made in 
Ref. [O] that local equilibration is achieved in 10.6 AGeV central Au+Au collisions at 
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t >10 fm/c within a volume of (5 fm) 3 . Apparently, deviations from angular isotropy can 
be explained by a nonspherical shape of the dense central region elongated in the beam 
direction (I7| . 

The interference between the bremsstrahlung photons from the two considered stages 
indicates to which extent it is justified to study the radiation from the different stages 
separately. From the comparison of Figs. [|a and |]c one can see that the relative contri- 
bution of this interference is of the order of 5 — 10%. The noticeable negative interference 
at small photon energies probably follows from a reacceleration of charges at late stages 
of the collision. It may serve as a signal of the final expansion of the compressed matter. 
Suppression of soft photons due to reacceleration of charges at late stages was found also 
in the Landau fluid-dynamical model |J. Due to the large interference the bremsstrahlung 
models considering only the stopping phase of a heavy-ion collision may underestimate 
the stopping power of nuclear matter if such models are used to fit the observed soft 
photon spectra. 

Figure [| represents the energy spectra of photons predicted with the SWM, DLM2 
and UrQMD model at the angle of maximum radiation given by Eq. ([]). One can see 
a qualitatively different behaviour of the spectra predicted by simple models and the 
UrQMD model at large u's. In the considered case of the AGS energy the SWM and DLM2 
noticeably underestimate the UrQMD predictions at uj > 150 MeV. Moreover, the UrQMD 
spectrum slightly increases with energy at uo > 200 MeV. Apparently such a behaviour also 
originates from the late stages of the reaction. 

As a result one may conclude that non-soft domains of the photon spectra can not 
be studied realistically without an explicit treatment of the late stages of a heavy-ion 
collision. On the other hand, the presence of the 7r°-background highly complicates the 
measurements of bremsstrahlung spectra at large photon energies. 

This is illustrated in Figure which shows the photon spectra predicted by the 
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URQMD model at different bombarding energies. In accordance with the general dis- 
cussion the photon spectra become more forward peaked with increasing initial energy. 
The background of electromagnetic decays of tt° mesons at the AGS energy exceeds the 
coherent photon yield already at uj > 50 MeV. At the SPS and RHIC energies the coherent 
bremsstrahlung remains important even for higher cj's (at small emission angles). 



V. CONCLUSIONS 

We have discussed general properties of the coherent bremsstrahlung spectra common 
to any model of heavy-ion collisions. The comparison of four different models shows that 
the photon spectra are sensitive to the dynamics of the projectile-target stopping only 
in the non-soft regime. We have shown, however, that even at higher photon energies 
different models may generate similar photon spectra. Photon bremsstrahlung produced 
in the stopping phase of a heavy-ion collision yields a dominant contribution to the soft 
photon spectra. As shown within the microscopic model the radiation of non-soft photons 
comes mainly from later (expansion) stages. The interference of this radiation with the 
emission from the initial stage is relatively small. 

The soft photon yield can provide some information about the degree of stopping. 
Soft photon spectra at ultrarelativistic bombarding energies, however, reveal a rather 
weak sensitivity to the final velocities of secondary particles. 

At the AGS and lower bombarding energies the background of 7r° decays is relatively 
large even for soft photons. According to our analysis, measurements of the soft photon 
spectra can hardly provide precise information about characteristics of nuclear stopping. 

It is shown by general analysis that oscillations of the energy spectra always appear 
in simple models assuming a finite time of the collision and, therefore, may be an artifact 
of a simplified description. On the other hand, using these oscillations to extract proper- 
ties of matter in nuclear collisions is only possible within certain simplified models [ 15] . 
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Moreover, these oscillations are strongly smoothened if one takes into account the photon 
radiation from the late expansion stages. 

The analysis of microscopic calculations performed within the UrQMD model shows 
that the shapes of the bremsstrahlung spectra as functions of the photon energy are 
strongly sensitive to the late stages of the reaction. In particular, this considerably com- 
plicates using spectral slopes as a tool to estimate the compression of matter in a heavy-ion 
collision. 
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FIG. 1. Trajectory representation of the shock wave model (left) and the degradation length 
model (right). Both calculations represent a central Au+Au collision at AGS energy (10.6 
AGeV). 
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FIG. 2. Angular distributions of soft photon bremsstrahlung (w = 1 MeV) predicted by the 
degradation length model with and without nuclear form factor. 
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FIG. 3. Angular distribution of bremsstrahlung spectra predicted by different models for a 
central 10.6 AGeV Au+Au collision at the photon energies oj = 5, 50, 150 and 250 MeV. The 
dots show the contribution of tt° — > 77 decays. 
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FIG. 4. Contributions to the bremsstrahlung spectra from different stages of a central 
Au+Au collision at AGS energy: a) deceleration stage (t < t cut ), b) expansion stage (t > t cut ), 
c) interference of a) and b). 
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FIG. 5. Energy spectra of bremsstrahlung photons in a 10.6 AGeV central Au+Au collision. 
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FIG. 6. Bremsstrahlung and tt° decay spectra predicted by the UrQMD model for central 
Au+Au collisions at the SIS, AGS, SPS and RHIC energies. 
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